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Abstract 



A new method for finding solutions of the nonlinear Shrodinger equation is 
proposed. Comutative multiplicative group of the nonlinear transformations, 
which operate on stationary localized solutions, enables a consideration of 
fractal subspaces in the solution space, stability and deterministic chaos. An 
increase of the transmission rate at the optical fiber communications can 
be based on new forms of localized stationary solutions, without significant 
change of input power. The estimated transmission rate is 50Gbit/s, for 
certain available soliton transmission systems. 



The propagation of pulsed light in an optical fiber can be described by 
the nonlinear Schrodinger equation, 

■ 9q(C,r) ld 2 q(£,r) 2 

where g(£, r) is a complex envelope function of the effective electric field 
amplitude and 

£ oc x, roc(t — x— ). (2) 

au; 

The higher order dispersion and the effect of fiber loss are neglected here |j] . 
We take 

q = q a e l ^y{r), (3) 
where y(r) is a real function, and get 

1 d 2 v o , , 

y - — -A - V = o. 4 

The solution of this equation 

yo(r) = (5) 
cosh g r 

describes the optical soliton. Its unchangeable shape is a property that makes 
it attractive for applaying to ultra high speed optical communications || |J . 
The equation (1) is completely integrable one. The inverse scattering trans- 
formation method || yields the general solutions of such nonlinear partial 
differential equations. Our aim is to propose here an alternative approach to 
the nonlinear Schrodinger equation and discus applicability of the obtained 
results to optical fiber communications. 

The equation (4) describes a stationary pulse in optical fiber. We take 
a localized solution y(r) of this equation and define the nonlinear operator 
H C1 : 

oo 

H c 1 y = Y, c jy^ ( 6 ) 
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where Cj are real coefficients. Does H ci y satisfy the equation (4)? The case 
y = yo is considered yet and the answer is positive 0. Putting H ci y into the 
equation (4), we find that H Cl y is actualy a solution of this equation if 

C2 3 = 0, (7) 
while C2j+i satisfy the recursive relation 

I 2j n-l 

c 2j+l = , ,A j( 2 3 ~ l)c 2j -l - C 2i+l-n Z! C «-fc c fc} ; (8) 

Z JU + -U n =2 fe=l 

where c\ is an arbitrary coefficient. Using the relations (6)- (8), with c\ = 1, 
we get 

= y- (9) 

In the following text, i? Cl will mean both the series (6) and the recursion (8) 
with (7). For a localized y(r) and a finite c±, convergence of the series (6) 
can be numerically tested. Our calculations yield that H Cl y is localized too. 
Therefore, using different values of ci, we are able to get uncountable many 
new localized solutions of the equation (4) from only one known localized 
solution (fig. 1). In the following text "the solution" will mean "the localized 
solution of the equation (4)". The solution value preciseness will be limited 
only by the number of calculated coeficients. The solution in form different 
from (6) does not exist. Each solution pair z(t) and y(r) must be in a relation 
z = H Cl y, with specific value of c\. 

Cl = lim 44- (10) 

r-»±oo y(j) 

Starting with a solution y{r) we can construct the complete solution space. 
There is an analogy to the superposition principle from linear theory. Ac- 
cording to the relation (10), a solution is determined by its asymptotics. 

The nonlinear Schrodinger equation has infinite number of symmetries 
corresponding to the conserved quantities: total energy, momentum, Hamil- 
tonian, ... JT[|. We find that there are actually uncountable conserved quan- 
tities. Let us consider the total energy only (for H Cl y): 



/oo 
( Cl y + c 3 y 3 + c 5 y 5 + ...)dr. (11) 
-oo 
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We can choose uncountable different values of C\ and use the relations (7) 
and (8). 



The relations (6)-(8) yield 



H ai H bl — H aibl . 



(12) 



Hence 



{# ci ;ci^0} 



(13) 



is the comutative multiplicative group of the nonlinear transformations (GNT). 
Group properties of the GNT originate from group properties of real numbers 
ci 7^ 0. For example, 



For definite coefficient c\ and solution y(r), we can construct a fractal sub- 
space of the solution space. The fractal subspace covers solutions of form 



In the phase plain, a fractal subspace is represented by a geometrical fractal 



For optical fiber communications it is important question whether small 
disturbations will destroy the information carrying pulses. Solution parame- 
ters, amplitude (pulse width) and velocity (frequency), are affected by various 
perturbations: outside produced noise, incoherence of the light source, fiber 
inhomogenities, absorption, amplifier noise, soliton interactions... It is the 
experimental fact that optical solitons (equation (5)) are unlikely to be de- 
stroyed by perturbations - they are very robust. We expect that at least the 
part of new solutions we have expresed here are actually stable. We are going 
to consider this problem theoretically, although it will be open until an ex- 
perimental verification. The GNT method enables the following statement: 
the stability of a solution y{r) is equivalent to the relation 




= #i. 



(14) 



H Cl H Cl ...H cl y. 



(15) 



(fig- 2). 



lim H 1+e y = y. 



(16) 



The relations (6)-(8) and (16) yield 



y(r) l< i. 



(17) 
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A localized solution of the equation (4) is stable one if and only if the relation 
(17) holds (fig. la,b). As well as for the KdV soliton 0, the classical ar- 
gument about the counterbalance between nonlinearity and dispersion is not 
sufficient to explain the stability. Consideration of the Lyapunov exponent, 

A(ci) = lim -In | |, (18) 

yields that deterministic chaos will appear at close packing of solitons, when 
Ci is large enough (fig. 3a). The deterministic chaos we can expect for 
ci > 2.4. Near c\ = 1, stability is exceptional (fig. 3b). 

New forms of localized stationary solutions of the nonlinear Schrodinger 
equation enable an increase of the transmission rate at the optical fiber com- 
munications, without significant change of input power. An information may 
be contained in the special form of soliton (fig. la,b). The known optical 
soliton, described by (5), is one of many possible stationary pulses. Let us 
consider an available soliton transmission system. If the fiber core cross sec- 
tional area is S — 60/im 2 , the carrier wavelength is A = 1.55/xra, the soliton 
pulse (equation (5)) width is r s = 25ps, the peak power is P m = 2.1mW, and 
the separation between two adjacent solitons is 3r s , then the transmission 
rate is 10Gbit/ s[8}. In the same transmission system, using stable pulses of 
form H ci y with different c\ (fig. la), the transmission rate will be greater. 
It becomes equal to 50Gbit/ s, at 40 photons resolution of energies. The new 
forms of stable solutions (fig. lb) make possible increase of the transmission 
rate in the same system. 

In conclusion, we have proposed the GNT method for solving of the non- 
linear Schrodinger equation. New forms of the stationary localized solutions, 
usable for an improvement of the optical fiber communications, are obtained. 

The authors would like to thank H.J.S.Dorren for useful discussion of the 
GNT method. This work was supported by the Soros Fund Open Society 
(Bosnia and Herzegovina) and the World University Service (Austria). 
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(a) Stable solutions with a = 0.2, 0.4, 0.6, 0.8, 1.0 

(b) Stable solutions with a = 1.4, 2.0, 2.39. 



Figure 1: (c)Unstable solutions with c\ = 2.41, 2.42. 
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Figure 2: Phase diagrams of the solutions (H Cl ) n y . (a)ci = 0.6, n — 1 to 8, 
(b)ci = 0.95, n = 1 to 5, (c)ci = 1.2, n = 1 to 5. 



(a)Lyapunov exponent 

Figure 3: (b) Exceptional stability near c\ = 1. 
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